Network analysis reveals a major role for 14q32 cluster miRNAs in determining transcriptional differences between IGHV-mutated and unmutated CLL.

Chronic lymphocytic leukaemia (CLL) cells can express unmutated (U-CLL) or mutated (M-CLL) immunoglobulin heavy chain (IGHV) genes with differing clinical behaviours, variable B cell receptor (BCR) signalling capacity and distinct transcriptional profiles. As it remains unclear how these differences reflect the tumour cells' innate pre/post germinal centre origin or their BCR signalling competence, we applied mRNA/miRNA sequencing to 38 CLL cases categorised into three subsets by IGHV mutational status and BCR signalling capacity. We identified 492 mRNAs and 38 miRNAs differentially expressed between U-CLL and M-CLL, but only 9 mRNAs and 0 miRNAs associated with BCR competence within M-CLL. Of the IGHV-associated miRNAs, (14/38 (37%)) derived from chr14q32 clusters where all miRNAs were co-expressed with the MEG3 lncRNA from a cancer associated imprinted locus. Integrative analysis of miRNA/mRNA data revealed pronounced regulatory potential for the 14q32 miRNAs, potentially accounting for up to 25% of the IGHV-related transcriptome signature. GAB1, a positive regulator of BCR signalling, was potentially regulated by five 14q32 miRNAs and we confirmed that two of these (miR-409-3p and miR-411-3p) significantly repressed activity of the GAB1 3'UTR. Our analysis demonstrates a potential key role of the 14q32 miRNA locus in the regulation of CLL-related gene regulation.

IL-33 Induces an Antiviral Signature in Mast Cells but Enhances Their Permissiveness for Human Rhinovirus Infection.

Mast cells (MCs) are classically associated with allergic asthma but their role in antiviral immunity is unclear. Human rhinoviruses (HRVs) are a major cause of asthma exacerbations and can infect and replicate within MCs. The primary site of HRV infection is the airway epithelium and MCs localise to this site with increasing asthma severity. The asthma susceptibility gene, IL-33, encodes an epithelial-derived cytokine released following HRV infection but its impact on MC antiviral responses has yet to be determined. In this study we investigated the global response of LAD2 MCs to IL-33 stimulation using RNA sequencing and identified genes involved in antiviral immunity. In spite of this, IL-33 treatment increased permissiveness of MCs to HRV16 infection which, from the RNA-Seq data, we attributed to upregulation of ICAM1. Flow cytometric analysis confirmed an IL-33-dependent increase in ICAM1 surface expression as well as LDLR, the receptors used by major and minor group HRVs for cellular entry. Neutralisation of ICAM1 reduced the IL-33-dependent enhancement in HRV16 replication and release in both LAD2 MCs and cord blood derived MCs. These findings demonstrate that although IL-33 induces an antiviral signature in MCs, it also upregulates the receptors for HRV entry to enhance infection. This highlights the potential for a gene-environment interaction involving IL33 and HRV in MCs to contribute to virus-induced asthma exacerbations.

Characterization of metabolic alterations of chronic lymphocytic leukemia in the lymph node microenvironment.

Altered metabolism is a hallmark of both cell division and cancer. Chronic lymphocytic leukemia (CLL) cells circulate between peripheral blood (PB) and lymph nodes (LNs), where they receive proliferative and prosurvival signals from surrounding cells. However, insight into the metabolism of LN CLL and how this may relate to therapeutic response is lacking. To obtain insight into CLL LN metabolism, we applied a 2-tiered strategy. First, we sampled PB from 8 patients at baseline and after 3-month ibrutinib (IBR) treatment, which forces egress of CLL cells from LNs. Second, we applied in vitro B-cell receptor (BCR) or CD40 stimulation to mimic the LN microenvironment and performed metabolomic and transcriptomic analyses. The combined analyses indicated prominent changes in purine, glucose, and glutamate metabolism occurring in the LNs. CD40 signaling mostly regulated amino acid metabolism, tricarboxylic acid cycle (TCA), and energy production. BCR signaling preferably engaged glucose and glycerol metabolism and several biosynthesis routes. Pathway analyses demonstrated opposite effects of in vitro stimulation vs IBR treatment. In agreement, the metabolic regulator MYC and its target genes were induced after BCR/CD40 stimulation and suppressed by IBR. Next, 13C fluxomics performed on CD40/BCR-stimulated cells confirmed a strong contribution of glutamine as fuel for the TCA cycle, whereas glucose was mainly converted into lactate and ribose-5-phosphate. Finally, inhibition of glutamine import with V9302 attenuated CD40/BCR-induced resistance to venetoclax. Together, these data provide insight into crucial metabolic changes driven by the CLL LN microenvironment. The prominent use of amino acids as fuel for the TCA cycle suggests new therapeutic vulnerabilities.

B-cell receptor dependent phagocytosis and presentation of particulate antigen by chronic lymphocytic leukemia cells.

Aim: T-helper cells could play an important role in the pathogenesis of chronic lymphocytic leukemia (CLL), a common B-cell neoplasm. Although CLL cells can present soluble antigens targeted from the B-cell receptor to T-helper cells via major histocompatibility complex (MHC) class II, antigens recognized by some CLL cells may be encountered in a particulate form. Here the ability of CLL cells to internalize and present anti-immunoglobulin M (IgM) beads as a model for the interaction of CLL cells with particulate antigens was investigated. Methods: The effect of anti-IgM beads on antigen presentation pathways was analyzed using RNA-seq and internalization of anti-IgM beads by primary CLL cells was investigated using confocal microscopy and flow cytometry. Antigen presentation was investigated by analyzing activation of a T-cell line expressing a T-cell receptor specific for a peptide derived from mouse κ light chains after incubating CLL cells with a mouse κ light chain-containing anti-IgM monoclonal antibody. Kinase inhibitors were used to characterize the pathways mediating internalization and antigen presentation. Results: Stimulation of surface IgM of CLL cells increased expression of the antigen presentation machinery and CLL cells were able to phagocytose anti-IgM beads. Internalization of anti-IgM beads was associated with MHC class II-restricted activation of cognate T-helper cells. Antigen presentation by CLL cells was dependent on activity of spleen tyrosine kinase (SYK) and phosphatidylinositol 3-kinase delta (PI3Kδ) but was unaffected by inhibitors of Bruton's tyrosine kinase (BTK). Conclusions: CLL cells can internalize and present antigen from anti-IgM beads. This capacity of CLL cells may be particularly important for recruitment of T-cell help in vivo in response to particulate antigens.

Insertion of atypical glycans into the tumor antigen-binding site identifies DLBCLs with distinct origin and behavior.

Glycosylation of the surface immunoglobulin (Ig) variable region is a remarkable follicular lymphoma-associated feature rarely seen in normal B cells. Here, we define a subset of diffuse large B-cell lymphomas (DLBCLs) that acquire N-glycosylation sites selectively in the Ig complementarity-determining regions (CDRs) of the antigen-binding sites. Mass spectrometry and X-ray crystallography demonstrate how the inserted glycans are stalled at oligomannose-type structures because they are buried in the CDR loops. Acquisition of sites occurs in ∼50% of germinal-center B-cell-like DLBCL (GCB-DLBCL), mainly of the genetic EZB subtype, irrespective of IGHV-D-J use. This markedly contrasts with the activated B-cell-like DLBCL Ig, which rarely has sites in the CDR and does not seem to acquire oligomannose-type structures. Acquisition of CDR-located acceptor sites associates with mutations of epigenetic regulators and BCL2 translocations, indicating an origin shared with follicular lymphoma. Within the EZB subtype, these sites are associated with more rapid disease progression and with significant gene set enrichment of the B-cell receptor, PI3K/AKT/MTORC1 pathway, glucose metabolism, and MYC signaling pathways, particularly in the fraction devoid of MYC translocations. The oligomannose-type glycans on the lymphoma cells interact with the candidate lectin dendritic cell-specific intercellular adhesion molecule 3 grabbing non-integrin (DC-SIGN), mediating low-level signals, and lectin-expressing cells form clusters with lymphoma cells. Both clustering and signaling are inhibited by antibodies specifically targeting the DC-SIGN carbohydrate recognition domain. Oligomannosylation of the tumor Ig is a posttranslational modification that readily identifies a distinct GCB-DLBCL category with more aggressive clinical behavior, and it could be a potential precise therapeutic target via antibody-mediated inhibition of the tumor Ig interaction with DC-SIGN-expressing M2-polarized macrophages.

DC-SIGN binding to mannosylated B-cell receptors in follicular lymphoma down-modulates receptor signaling capacity.

In follicular lymphoma (FL), surface immunoglobulin (sIg) carries mandatory N-glycosylation sites in the variable regions, inserted during somatic hypermutation. These glycosylation sites are tumor-specific, indicating a critical function in FL. Added glycan unexpectedly terminates at high mannose (Mann) and confers capability for sIg-mediated interaction with local macrophage-expressed DC-SIGN lectin resulting in low-level activation of upstream B-cell receptor signaling responses. Here we show that despite being of low-level, DC-SIGN induces a similar downstream transcriptional response to anti-IgM in primary FL cells, characterized by activation of pathways associated with B-cell survival, proliferation and cell-cell communication. Lectin binding was also able to engage post-transcriptional receptor cross-talk pathways since, like anti-IgM, DC-SIGN down-modulated cell surface expression of CXCR4. Importantly, pre-exposure of a FL-derived cell line expressing sIgM-Mann or primary FL cells to DC-SIGN, which does not block anti-IgM binding, reversibly paralyzed the subsequent Ca2+ response to anti-IgM. These novel findings indicate that modulation of sIg function occurs in FL via lectin binding to acquired mannoses. The B-cell receptor alternative engagement described here provides two advantages to lymphoma cells: (i) activation of signaling, which, albeit of low-level, is sufficient to trigger canonical lymphoma-promoting responses, and (ii) protection from exogenous antigen by paralyzing anti-IgM-induced signaling. Blockade of this alternative engagement could offer a new therapeutic strategy.

Clinical significance of TP53, BIRC3, ATM and MAPK-ERK genes in chronic lymphocytic leukaemia: data from the randomised UK LRF CLL4 trial.

Despite advances in chronic lymphocytic leukaemia (CLL) treatment, globally chemotherapy remains a central treatment modality, with chemotherapy trials representing an invaluable resource to explore disease-related/genetic features contributing to long-term outcomes. In 499 LRF CLL4 cases, a trial with >12 years follow-up, we employed targeted resequencing of 22 genes, identifying 623 mutations. After background mutation rate correction, 11/22 genes were recurrently mutated at frequencies between 3.6% (NFKBIE) and 24% (SF3B1). Mutations beyond Sanger resolution (<12% VAF) were observed in all genes, with KRAS mutations principally composed of these low VAF variants. Firstly, employing orthogonal approaches to confirm <12% VAF TP53 mutations, we assessed the clinical impact of TP53 clonal architecture. Whilst ≥ 12% VAF TP53mut cases were associated with reduced PFS and OS, we could not demonstrate a difference between <12% VAF TP53 mutations and either wild type or ≥12% VAF TP53mut cases. Secondly, we identified biallelic BIRC3 lesions (mutation and deletion) as an independent marker of inferior PFS and OS. Finally, we observed that mutated MAPK-ERK genes were independent markers of poor OS in multivariate survival analysis. In conclusion, our study supports using targeted resequencing of expanded gene panels to elucidate the prognostic impact of gene mutations.

Clinical significance of DNA methylation in chronic lymphocytic leukemia patients: results from 3 UK clinical trials.

Chronic lymphocytic leukemia patients with mutated immunoglobulin heavy-chain genes (IGHV-M), particularly those lacking poor-risk genomic lesions, often respond well to chemoimmunotherapy (CIT). DNA methylation profiling can subdivide early-stage patients into naive B-cell-like CLL (n-CLL), memory B-cell-like CLL (m-CLL), and intermediate CLL (i-CLL), with differing times to first treatment and overall survival. However, whether DNA methylation can identify patients destined to respond favorably to CIT has not been ascertained. We classified treatment-naive patients (n = 605) from 3 UK chemo and CIT clinical trials into the 3 epigenetic subgroups, using pyrosequencing and microarray analysis, and performed expansive survival analysis. The n-CLL, i-CLL, and m-CLL signatures were found in 80% (n = 245/305), 17% (53/305), and 2% (7/305) of IGHV-unmutated (IGHV-U) cases, respectively, and in 9%, (19/216), 50% (108/216), and 41% (89/216) of IGHV-M cases, respectively. Multivariate Cox proportional analysis identified m-CLL as an independent prognostic factor for overall survival (hazard ratio [HR], 0.46; 95% confidence interval [CI], 0.24-0.87; P = .018) in CLL4, and for progression-free survival (HR, 0.25; 95% CI, 0.10-0.57; P = .002) in ARCTIC and ADMIRE patients. The analysis of epigenetic subgroups in patients entered into 3 first-line UK CLL trials identifies m-CLL as an independent marker of prolonged survival and may aid in the identification of patients destined to demonstrate prolonged survival after CIT.

Characterization of Somatically-Acquired Copy Number Alterations in Chronic Lymphocytic Leukaemia Using Shallow Whole Genome Sequencing.

Shallow whole genome sequencing (sWGS) is a simple, robust, and cost-effective technique recently optimized for the identification of copy number aberrations (CNAs) in tumor samples. This multiplexed methodology sequences 50 bp from one end of the DNA molecule, generating ˜0.1× coverage, and utilizes the observed sequence depth across the genome to infer copy number. It is amenable to low quantities of input DNA, sequencing costs are modest, processing is compatible with low-output instruments, and downstream analysis is simplified by the use of freely available bioinformatics tools and a data analysis package written especially for the analysis of sWGS data. It is the aim of this chapter to introduce the fundamental concepts of sWGS and to provide an overview of the steps involved in a successful sWGS experiment.

Human Papillomavirus DNA Methylation Predicts Response to Treatment Using Cidofovir and Imiquimod in Vulval Intraepithelial Neoplasia 3.

Purpose: Response rates to treatment of vulval intraepithelial neoplasia (VIN) with imiquimod and cidofovir are approximately 57% and 61%, respectively. Treatment is associated with significant side effects and, if ineffective, risk of malignant progression. Treatment response is not predicted by clinical factors. Identification of a biomarker that could predict response is an attractive prospect. This work investigated HPV DNA methylation as a potential predictive biomarker in this setting.Experimental Design: DNA from 167 cases of VIN 3 from the RT3 VIN clinical trial was assessed. HPV-positive cases were identified using Greiner PapilloCheck and HPV 16 type-specific PCR. HPV DNA methylation status was assessed in three viral regions: E2, L1/L2, and the promoter, using pyrosequencing.Results: Methylation of the HPV E2 region was associated with response to treatment. For cidofovir (n = 30), median E2 methylation was significantly higher in patients who responded (P ≤ 0.0001); E2 methylation >4% predicted response with 88.2% sensitivity and 84.6% specificity. For imiquimod (n = 33), median E2 methylation was lower in patients who responded to treatment (P = 0.03; not significant after Bonferroni correction); E2 methylation <4% predicted response with 70.6% sensitivity and 62.5% specificity.Conclusions: These data indicate that cidofovir and imiquimod may be effective in two biologically defined groups. HPV E2 DNA methylation demonstrated potential as a predictive biomarker for the treatment of VIN with cidofovir and may warrant investigation in a biomarker-guided clinical trial. Clin Cancer Res; 23(18); 5460-8. ©2017 AACR.

Exome Sequencing in Classic Hairy Cell Leukaemia Reveals Widespread Variation in Acquired Somatic Mutations between Individual Tumours Apart from the Signature BRAF V(600)E Lesion.

In classic Hairy cell leukaemia (HCLc), a single case has thus far been interrogated by whole exome sequencing (WES) in a treatment naive patient, in which BRAF V(600)E was identified as an acquired somatic mutation and confirmed as occurring near-universally in this form of disease by conventional PCR-based cohort screens. It left open however the question whether other genome-wide mutations may also commonly occur at high frequency in presentation HCLc disease. To address this, we have carried out WES of 5 such typical HCLc cases, using highly purified splenic tumour cells paired with autologous T cells for germline. Apart from BRAF V(600)E, no other recurrent somatic mutation was identified in these HCLc exomes, thereby excluding additional acquired mutations as also prevalent at a near-universal frequency in this form of the disease. These data then place mutant BRAF at the centre of the neoplastic drive in HCLc. A comparison of our exome data with emerging genetic findings in HCL indicates that additional somatic mutations may however occur recurrently in smaller subsets of disease. As mutant BRAF alone is insufficient to drive malignant transformation in other histological cancers, it suggests that individual tumours utilise largely differing patterns of genetic somatic mutations to coalesce with BRAF V(600)E to drive pathogenesis of malignant HCLc disease.

Massive parallel IGHV gene sequencing reveals a germinal center pathway in origins of human multiple myeloma.

Human multiple myeloma (MM) is characterized by accumulation of malignant terminally differentiated plasma cells (PCs) in the bone marrow (BM), raising the question when during maturation neoplastic transformation begins. Immunoglobulin IGHV genes carry imprints of clonal tumor history, delineating somatic hypermutation (SHM) events that generally occur in the germinal center (GC). Here, we examine MM-derived IGHV genes using massive parallel deep sequencing, comparing them with profiles in normal BM PCs. In 4/4 presentation IgG MM, monoclonal tumor-derived IGHV sequences revealed significant evidence for intraclonal variation (ICV) in mutation patterns. IGHV sequences of 2/2 normal PC IgG populations revealed dominant oligoclonal expansions, each expansion also displaying mutational ICV. Clonal expansions in MM and in normal BM PCs reveal common IGHV features. In such MM, the data fit a model of tumor origins in which neoplastic transformation is initiated in a GC B-cell committed to terminal differentiation but still targeted by on-going SHM. Strikingly, the data parallel IGHV clonal sequences in some monoclonal gammopathy of undetermined significance (MGUS) known to display on-going SHM imprints. Since MGUS generally precedes MM, these data suggest origins of MGUS and MM with IGHV gene mutational ICV from the same GC B-cell, arising via a distinctive pathway.

Human papillomavirus type 16 L1/L2 DNA methylation shows weak association with cervical disease grade in young women.

BACKGROUND: Persistent infection with human papillomavirus (HPV) type 16 causes the majority of cervical cancers. Genital HPV infection is very common, but neoplastic progression is uncommon. There is an urgent need for biomarkers associated with cervical neoplasia, to enable triage of women who test positive for HPV. OBJECTIVES: To assess the ability of quantitative measurement of HPV16 DNA methylation to separate samples of different cytological and histological grades from young women, among whom rates of HPV infection are high. STUDY DESIGN: DNA methylation was quantified by pyrosequencing of bisulphite converted DNA from liquid based cytology samples from 234 women (mean age 20.6 years) who tested positive for HPV16 and showed varying degrees of neoplasia. Methylation was assessed at CpGs in the HPV E2 and L1/L2 regions. RESULTS: The performance of methylation-based classifiers was assessed by ROC curve analyses. The best combination of CpGs (5600 and 5609) achieved AUCs of 0.656 (95% CI=0.520-0.792) for separation of cytologically normal and severely dyskaryotic samples, and 0.639 (95% CI=0.547-0.731) for separation of samples with or without high-grade neoplasia (CIN2+/-). CONCLUSIONS: The data are consistent with HPV L1/L2 methylation being a marker of the duration of infection in a specific host. Assessment of HPV DNA methylation is hence a promising biomarker to triage HPV-positive cytology samples, but may have limited utility in young women. Future studies assessing the likely utility of HPV DNA methylation as a potential triage biomarker must take account of women's age.

Increased methylation of Human Papillomavirus type 16 DNA correlates with viral integration in Vulval Intraepithelial Neoplasia.

BACKGROUND: Methylation of HPV16 DNA is a promising biomarker for triage of HPV positive cervical screening samples but the biological basis for the association between HPV-associated neoplasia and increased methylation is unclear. OBJECTIVES: To determine whether HPV16 DNA methylation was associated with viral integration, and investigate the relationships between viral DNA methylation, integration and gene expression. STUDY DESIGN: HPV16 DNA methylation, integration and gene expression were assessed using pyrosequencing, ligation-mediated PCR and QPCR, in biopsies from 25 patients attending a specialist vulval neoplasia clinic and in short-term clonal cell lines derived from vulval and vaginal neoplasia. RESULTS: Increased methylation of the HPV16 L1/L2 and E2 regions was associated with integration of viral DNA into the host genome. This relationship was observed both in vivo and in vitro. Increased methylation of E2 binding sites did not appear to be associated with greater expression of viral early genes. Expression of HPV E6 and E7 did not correlate with either integration state or increased L1/L2 methylation. CONCLUSIONS: The data suggest that increased HPV DNA methylation may be partly attributable to viral integration, and provide a biological rationale for quantification of L1/L2 methylation in triage of HPV positive cervical screening samples.

mRNA sequencing of novel cell lines from human papillomavirus type-16 related vulval intraepithelial neoplasia: consequences of expression of HPV16 E4 and E5.

Vulval intraepithelial neoplasia is a precursor of vulval cancer and is commonly caused by infection with Human Papillomavirus (HPV). Development of topical treatments for vulval intraepithelial neoplasia requires appropriate in vitro models. This study evaluated the feasibility of primary culture of vulval intraepithelial neoplasia biopsy tissue to produce cell lines for use as in vitro models. A potentially immortal cell line was produced which gave rise to three monoclonal lines. These lines were characterized for HPV genomic integration and for viral gene expression using ligation-mediated PCR and quantitative PCR. Distinct patterns of viral integration and gene expression were observed among the three lines. Integration and expression data were validated using deep sequencing of mRNA. Gene ontology analyses of these data also demonstrated that expression of the HPV16 E4 and E5 proteins resulted in substantial changes in the composition of the cell membrane and extracellular space, associated with alterations in cell adhesion and differentiation. These data illustrate the diverse patterns of HPV gene expression potentially present within a single lesion. The derived cell lines provide useful models to investigate the biology of vulval intraepithelial neoplasia and the interactions between different HPV gene products and potential therapeutic agents.

Quantitative measurement of Human Papillomavirus type 16 L1/L2 DNA methylation correlates with cervical disease grade.

BACKGROUND: Persistent infection with Human Papillomavirus (HPV) type 16 causes the majority of cervical cancers. Genital HPV infection is very common, but neoplastic progression is uncommon. There is an urgent need to identify biomarkers associated with cervical neoplasia that can be used to triage women who test positive for HPV. OBJECTIVES: To assess the ability of quantitative measurement of HPV16 DNA methylation to separate samples of different cytology grades and cervical cancers, and determine which of the assessed regions of the HPV genome and individual CpGs are most informative. STUDY DESIGN: DNA methylation was quantified by pyrosequencing of bisulphite converted DNA from liquid based cytology samples from 17 women with normal cytology and 20 women with severe dyskaryosis, and from fixed tissue from 24 women with cervical cancer. Methylation was assessed in the HPV Long Control Region (LCR), E2 and L1/L2 regions. RESULTS: In cervical cancers, increased HPV DNA methylation was present in all regions. Increased methylation was also observed in severely dyskaryotic relative to normal samples, but only in the E2 and L1/L2 regions. The ability of methylation based classifiers to separate the three classes of material was assessed by ROC curve analyses. The best separation between normal and dyskaryotic samples was achieved by assessment of the L1/L2 CpGs at nucleotide positions 5600 and 5609 (AUC=0.900, 95% CI: 0.793-1). CONCLUSIONS: This study demonstrates the potential of quantification of HPV DNA methylation as a biomarker of cervical neoplasia. An algorithm considering methylation at specific L1/L2 CpGs appeared the most promising model.

Human papillomavirus type distribution in vulval intraepithelial neoplasia determined using PapilloCheck DNA Microarray.

Vulval intraepithelial neoplasia is a precursor of vulval carcinoma, and is frequently associated with human papillomavirus (HPV) infection. Estimates of HPV prevalence in vulval intraepithelial neoplasia vary widely in the UK. The objective of this study was to assess HPV infection in a sample of women with vulval intraepithelial neoplasia, confirmed histologically, and determine the proportion of disease associated with HPV types targeted by prophylactic HPV vaccines. HPV infection was assessed in biopsies from 59 patients using the Greiner Bio-One PapilloCheck® DNA chip assay. Valid results were obtained for 54 cases. HPV infection was present in 43 of the 54 cases (79.6%: 95% CI 67.1-88.2%). The most common HPV types were HPV 16 (33/54: 61.1%), HPV 33 (8/54: 14.8%), HPV 6 (5/54: 9.3%), and HPV 42 (3/54: 5.6%). The mean age of HPV positive women was significantly less than the mean age of HPV negative women. This is the largest UK series of vulval intraepithelial neoplasia in which HPV type has been investigated, and 34/54 (63.0%, 95% CI: 49.6-78.6%) cases were associated with HPV 16/18, which are targeted by current prophylactic HPV vaccines.